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ABSTRACT 



Context. Central active galactic nuclei (AGN) are supposed to play a key role in the evolution of their host galaxies. In particular, the 
dynamical and physical properties of the gas core must be affected by the injected energy. 

Aims. Our aim is to study the effects of an AGN on the dark matter profile and on the central stellar light distribution in massive early 
type galaxies. 

Methods. By performing self-consistent N-body simulations, we assume in our analysis that periodic bipolar outbursts from a central 
AGN can induce harmonic oscillatory motions on both sides of the gas core. 

Results. Using realistic AGN properties, we find that the motions of the gas core, driven by such feedback processes, can flatten the 
dark matter and/or stellar profiles after 4-5 Gyr. These results are consistent with recent observational studies that suggest that most 
giant elliptical galaxies have cores or are "missing light" in their inner part. Since stars behave as a "collisionless" fluid similar to 
dark matter, the density profile both of stars and dark matter should be affected in a similar way, leading to an effective reduction in 
the central brightness. 

Key words, methods: N-body simulations - galaxies: structure - galaxies: active 



1. Introduction 



The cold dark matter (CDM) paradigm (Cole et al. l2005l and ref- 
erences therein) has led to a successful explanation of the large- 
scale structure in the galaxy distribution on scales 0.02 < k i 
0. 15h Mpc . The CDM power spectrum on these scales derived 
from large redshift surveys such as, for instance, the Anglo- 
Australian 2-degree Field Galaxy Redshift Survey (2dFGRS), 
is also consistent with the Lyman-o' forest data in the redshift 
range 2< z <4 (Croft et al. [2tjU2b . 

In spite of these impressive successes, there are still dis- 
crepancies between simulations and observations on scales i 
1 Mpc. We may mention the large number of sub-L» subhalos 
present in simulations but not observed (Kauffmann, White & 
Guiderdoni [T993l Moo re et al. |1999al Klypin et al. [19"9"9l but see 
Belokurov et al. 120061 for an observational update), and the ex- 
cess of massive early-type galaxies undergoing "top-down" as- 
sembly with high inferred specific star formation rates relative to 
predictions of the hierarchical scenario (Glazebrook et al. 2004; 
Cimatti, Daddi & Renzini 2006). Here we examine the sharp 
central density cusp predicted by simulations in dark matter ha- 
los and not seen in the rotation curves of low surface brightness 
galaxies (de Blok et al.|2p01 ) or in bright spiral galaxies (Palunas 
& Williams|2000l Salucci & Burkert|2000l Gentile [20041). Little 
attention has been given to the corresponding situation in mas- 
sive early-type galaxies, the focus of the present analysis. 

Our model is presented as a toy model to justify an in depth 
study. It is in the same spirit as that developed by Mashchenko, 
Couchman & Wadsley (2006]) to account for cores in dwarf 
galaxies. Indeed it has a similar problem: where has the gas 



gone? But we stress that there are many ways of accounting for 
this and other objections (we provide more details below). The 
point we wish to make is that central cores in massive ellipti- 
cals are possible tracers of violent AGN activity, and provide 
an alternative to binary black hole scouring. If AGN are power- 
ful sources of negative feedback that cleans out the gas supply 
to avoid forming excessively luminous and blue ellipticals, our 
model is a plausible and even inevitable precursor. These cores 
are seen in the stellar component, and are likely present in the 
underlying (although subdominant) dark matter distribution. 

Navarro, Frenk & White d 19961 119971 ) have shown that the 
spherically-averaged density profiles of simulated halos can be 
fitted by a simple analytical function, depending on a character- 
istic density p* and a characteristic radius r s : 
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This profile (dubbed the NFW-profile) is steeper than that of an 
isothermal sphere at large radii and shallower close to the center. 
Steeper profiles in the central regions (a ~ - 1 .5) have been found 
in high resolution simulations (Moore et al. 1999b; Ghigna et al. 
120001 Fukushige & Makino |20TJTl i. 

While the "universality" of density profiles of dark halos is 
still a matter of debate, possibly depending on the merger history 
(Klypin et al. 120011 Ricotti EU03l Boylan- Kolchin & Ma l2"004i l 
or on the initial conditions (Ascasibar et al. 2004), it could nev- 
ertheless be an important key to understanding the mechanism(s) 
by which these systems relax and attain equilibrium. In a colli- 
sionless self-gravitating system, collective mechanisms such as 
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violent relaxation and/or phase mixing certainly play a major 
role in the relaxation of halos, but density profiles may also be 
affected by gravitational scattering of dark matter particles in 
substructures present inside halos (Ma & Bertschinger 2004a). 
In fact, recent numerical experiments indicate that this diffusion 
process may alter significantly the inner density profile of ha- 
los, producing a flattening of the original profile within a few 
dynamical time scales (Ma & Boylan-Kolchin 2004b). The most 
detailed study (Ricotti 2003 ) of N-body simulations at different 
scales concluded that galaxies have shallower density profiles 
than clusters. 

Central cusps seen in simulations could also be understood 
as a consequence of the inflow of low-entropy material and 
therefore contain information on the relic entropy of dark mat- 
ter particles. In fact, high resolution simulations of galaxy-sized 
CDM halos indicate an increase of the coarse-grained phase- 
space density Q (defined as the ratio between the density and 
the cube of 1-D velocity dispersion in a given volume, e.g., 
Q = p/cr 3 ) towards the center (Taylor & Navarro 200 1)- Similar 
results were obtained by Rasia, Tormen & Moscardini (2003), 
who also obtained a power-law variation, e.g., Q oc for 
cluster-size halos, with f3 quite close to the value found by Taylor 
& Navarro (USD, namely, /? * 1.87. 

Several solutions have been proposed to explain such dis- 
crepancies between observations and numerical simulations. For 
example, the dark matter (DM) can be "heated" by the baryons 
by dynamical friction due to self-gravitating gas clouds orbiting 
near the center of the galaxy (El-Zant, Shlosman & Hoffman 
|200T| El-Zant et al. |20T)4l ), by the evolution of a stellar bar 
(Weinberg & Katz 2002, Holley-Bockelmann, Weinberg & Katz 
12005 1 Sellwood 2006 ), by the radiation recoil by a black hole 
(Merritt et al. 2004) or by random bulk gas motions driven 
by supernovae feedback, recently suggested by Mashchenko, 
Couchman & Wadsley (2006). Other mechanisms have been 
proposed such as the the transfer of angular momentum from 
baryonic to dark matter (Tonini, Lapi & Salucci2006 ) or the ex- 
pulsion of a large fraction of the gas due to feedback activities, 
causing the dark matter to expand (Gnedin & Zhao 120021 . 

Most of the previous studies consider the scales of dwarf 
galaxies or those of galaxy clusters in order to interpret the ob- 
served DM profiles in these objects which are thought to be dark 
matter dominated in their inner regions. In the present paper, 
we focus our attention on the scale of a typical giant elliptical 
galaxy, and we investigate the effects of a central AGN in its 
core. Similar but more extreme feedback ideas have been pro- 
posed for disk galaxies such as the Milky Way, where massive 
early outflows were invoked that resulted in homogenizing and 
heating the inner dark halo (Binney, Gerhard & Silk l200fl ). 

Black holes and associated AGN feedback are generally con- 
sidered to be active in the centers of giant galaxies. They must 
represent huge sources of energy and therefore are thought to 
play a key role in the formation of bright ellipticals (Silk 1998 ). 
They can also be responsible for heating the gas core to reg- 
ulate the cooling flows in both elliptical galaxies (Best et al. 
2006a) and brightest group and cluster galaxies (Dunn & Fabian 
2006; Best et al. 2006b). Although over the past few years, im- 
portant progress have been done with high spectral resolution 
using HST, Fuse, Chandra and XMM-Newton, the nature of the 
interaction between the AGNs and the surrounding intergalactic 
medium is still poorly understood. However, observations from 
UV and X-ray absorption lines of bright AGNs seem to suggest 
their outbursts can induce sonic motions to gas clouds and lead 
to important biconical mass outflows (see for instance Hutchings 
et al. [T9981 Bridges & Irwin [T9981 Churazov et al. [20041 and 



references therein). In the present work, we assume that each 
AGN outburst can induce a bulk motion to the gas core by en- 
ergy transfer. We then explore the efficiency of this mechanism 
of gas core bulk motions for flattening the central DM and stel- 
lar cusp. This paper is organized as follows: in section 2, we 
describe our toy model, we present the results in section 3, and 
in section 4, our main conclusions are summarized. 

2. Toy model 

We first develop a toy model which consists of building an N- 
body realization of an isolated, equilibrium model galaxy. A di- 
rect method uses the distribution function (DF), which provides 
the relative probability of a star to have a certain position and ve- 
locity, and then use a Monte-Carlo sampling of the DF to gener- 
ate the N-body realization. Thus, assuming spherical symmetry, 
we generate a spherical equilibrium DM halo with a Hernquist 
density profile (Hernquist |1990l ): 



where M^ m is the total mass of dark matter and a the scale ra- 
dius. In our fiducial model, we use Md„, = 10 13 h~ l M Q and a = 
77.0 h~ l kpc. For comparison, this profile is identical in the in- 
ner parts to an NFW-profile with virial radius r2oo = 445 h~ x kpc 
(r2oo defines the sphere within which the mean density is equal 
to 200 times the critical density) and scale radius r s = 44.5 h 
kpc respectively (see for instance relation 2 between a and r s in 
Springel, Di Matteo & Hernquist 2005a). These 2 profiles are 
compared in figure Q] The choice of using an Hernquist profile 
is motivated by two main reasons. Firstly, the density declines 
faster than an NFW-profile in the outer parts of the halo. Thus, 
for an isolated halo, truncation at the virial radius is not needed. 
Furthermore, the Hernquist profile has an analytical expression 
for both the distribution function and the velocity dispersion 
(Hernquist 119901 ), which is not the case for an NFW-profile. 
By proceeding in this way, we try to minimize undesirable ef- 
fects which can lead to unstable equilibrium. Conventionally, 
the scale radius is given in terms of the concentration parame- 
ter c = r v /r s . For our fiducial model, we use c - 10, in good 
agreement with values found in previous cosmological N-body 
simulations (Bullock et al. 12001 at Dolag et al. |2004| ). 

The gas component follows a Hernquist profile as well with 
a total gas mass M gas constrained by the baryonic fraction fi, = 
M gas /{Mdm + Mg as ). Its initial temperature profile is chosen so 
that the gas is initially in hydrostatic equilibrium. 

To model the effects of the central AGN on the gas compo- 
nent, we use a simple representation where the outflow takes 
the form of a collimated jet which causes the gas to expand 
by compression. Since the resulting detailed shape of the gas 
is complicated and poorly known, we use the following simpli- 
fied model of spherical symmetry (spherical symmetry makes 
the model well-defined with the minimum number of free pa- 
rameters). We first define a gas core radius which limits the re- 
gion affected by the AGN activity. We assume its value to be 
30h~ l kpc corresponding to less than 10% of the virial radius 
and includes a few percent of the total mass of the system. Then 
the gas core is split into 2 parts and due to periodic AGN out- 
bursts, each part has an independent harmonical motion on op- 
posite sides of the AGN. The frequency and the amplitude of the 
oscillations depend on the frequency of the AGN activity, on the 
energy transfer between the AGN and the gas, and on the phys- 
ical properties of the gas (viscosity, temperature, local density, 
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Fig. 1. DM density profile of our unperturbed model. The dots 
are estimates of the density at 5h~ l Gyr while the solid line 
is the expected Hernquist profile whose relevant parameters are 
discussed in the text. For comparison, the "equivalent" NFW- 
profile is also shown (dashed line). 



etc.). To simplify matters, we assume that i) the frequency of 
the AGN emission matches the frequency of the oscillations of 
the gas and ii) the amplitude of motions of a specific model cor- 
responds to a fractional displacement of the scale radius. Finally, 
we have also studied, in a second experiment, how our results are 
affected when a stellar component is included. 



3. Results 

3.1. Adiabatic gas 

We carry out each simulation with the public code GADGET2 
(Springel 2005b) where 5 x 10 5 dark matter particles and 5 x 10 5 
gas particles have been used. The value of the softening length is 
1 h kpc for all experiments and the baryonic fraction is taken 
to fb = 0.13. Although the gas is supposed to have a central 
core, our representation of a Hernquist profile will not affect our 
results since the relevant parameter is the total mass of the core 
gas. In our study, its resulting value is taken to be 1 .2 x 10 U M Q , 
i.e. 1% of the total mass of the system. 

The strongest effects are expected to occur when the veloc- 
ities of the gas particles are close to the velocity dispersion of 
dark matter. Then, in our fiducial model, we consider that each 
injection of energy by the AGN produces a gas bulk motion with 
an initial velocity of V gas = 260 km/s, corresponding to the ve- 
locity dispersion of dark matter at radius 10 - 20 hr x kpc. Since 
the escape velocity at r — 30 h~ x kpc is about 900 km/s, this en- 
sures that the two moving parts of the gas core are always bound 
to the system. Then, we let the system evolve during 5 h~ x Gyr 
which corresponds to about 20 times the dynamical time of the 
gas core. Figure Q] shows the DM density profile at t = 5h~ l Gyr 
for an unperturbed system (Vgas = 0). The whole system seems 
to be in a very stable equilibrium since the estimated densities at 
different radii match well those of the initial density profile. 

To study the time evolution of the degree of cuspiness of 
each model, we use a simple estimator defined by 



A = 

C = 0.00161 + 0.048 
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Fig. 2. Time evolution of the parameter <x for the model with 
Vg as = 260 km/s and for 6 different amplitudes A. For each ex- 
periment, the solid line is the best fit to the data. 



where p e xp(n) is the estimated density at the radius r, < 50 
kpc compared to the theoretical value p t heo(rd. N represents the 
number of measures used for the estimation of cr. Typically, we 
took N - 50 points separated within logarithmic bins of equal 
width. It is worth mentioning than for each time-step we have es- 
timated the mass center of the DM component by using a friend- 
of-friends algorithm (Davis et al. |1985l >. The value of the link- 
ing length is 0.05 in units of the mean interparticle separation, 
corresponding to evaluating the center of mass of the 5000 DM 
particles. Figure [2] shows the evolution of our estimator cr for 6 
specific amplitudes A of the oscillations. First, when the system 
is at rest, namely V = and A = 0, cr is constant and close to 
as expected. This experiment is used to check the stability of 
the equilibrium of our system and to estimate the maximum nu- 
merical error. The positive values are mostly due to the lack of 
resolution in the very inner parts of the halo (r, < 0.03 r s ) which 
can lead to underestimates or overestimates of the density. For 
both A = 1/10 r s and A = 1/5 r s , in spite of erratic fluctuations, 
on average, a small increase of cr is observed. On the contrary, 
for A = 1/2 r s , 3/4 r s and r s , cr increases faster but no significant 
differences can been seen between the 3 scenarios. 

However, two important remarks can be made. Firstly, the 
strongest effects occur when A = 1 /2r s where the DM particles 
seem to enter into resonance. Secondly, the value of cr oscil- 
lates with the same frequency as the gas oscillations. This phe- 
nomenon can be clearly seen for A = r s where the period of the 
oscillations is 0.52 ft -1 Gyr. This suggests that the shape of the 
inner profile continuously evolves over time through different 
phases. 
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Fig. 3. DM density profile for V gas - 260 km/s and A = 1 /2 r s 
plotted at three different times, t — 1 (dashed line), t — 2 (dotted 
line) and t - Ah~ x Gyr. The solid line represents the expected 
profile in the absence of any perturbation. 
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Fig. 4. Time evolution of the parameter cr for A = 1/2 r s and for 
Vg as = 50 km/s (right panel) or V gas = 500 km/s (right panel). 



Next, we consider the case A — l/2r s . Figure [3] represents 
the evolution of the DM density profile at t — 1,2 and 4 hr x Gyr. 
We notice that the DM cusp is progressively flatten. At t — 4 h' 1 
Gyr, the DM profile can be decomposed into a flat central part 
(r < 0.lr s ~ 5h~ l kpc) and a steeper outer part until it converges 
to the original slope at r ~ 0.7r s . In this case, the derived density 
profile cannot be well fitted by a Burkert profile (Burkert [1995l l. 

Figure|4]shows the evolution of the time evolution of <x when 
the velocity of the gas is either small (V gas = 50 km/s) or high 
(Vgas = 500 km/s) compared to the velocity dispersion of the 
dark matter. In both case, an increase of cr is observed but the 
effects are quite small, in good agreement with what we expect. 
It is interesting to notice that for V gas = 50 km/s, the three peaks 
correspond to 1, 2 and 3 oscillation periods respectively (the pe- 
riod of the oscillation is about 1.4 h~ l Gyr). At these specific 
times, the effects on the DM profile are negligible. 

3.2. Including a stellar distribution 

One way to improve our model is to include the formation of 
the elliptical galaxy itself; dissipation of the gas (from radia- 
tive cooling processes) and subsequent star formation leads to 
a steeper initial DM density profile due to adiabatic contraction. 
We achieve this by building new initial conditions resulting from 
the merger of 4 sub-haloes; this is done in order to boost the star 
formation rate (SFR) and form a giant galaxy in a relatively short 
period. 



Each sub-halo has a mass of 2.5 x 1O 12 M and the same 
density profiles as discussed in section 2. Additionally, we give 
each sub-halo angular momentum by setting the spin parameter, 
A = J\E\ l l 2 IGM 5 ! 2 » 0.1 (where J is the angular momentum, E 
is the total energy of the halo and M is its mass), and the spe- 
cific angular momentum distribution, j(r) oc r, as expected from 
cosmological A^-body simulations (Bullock et al. 2001b). The 4 
respective centers of mass are initially placed at rest in a cross 
configuration, of length 400 hr x kpc. Each opposing pair has op- 
positely orientated angular momentum vectors in order to avoid 
an excessive amount of rotation after the merger event. 

Before repeating our previous experiment, we first let the 
above configuration evolve for 6.5 h~ l Gyr. During the first 1.5 
h~ l Gyr (which captures the merger event), we turn on radia- 
tive cooling and star formation, then allow the system to relax 
for the remaining time. Cooling and star formation were intro- 
duced in the manner of Katz, Weinberg & Hernquist 1996 The 
cooling prescription allowed each gas particle with T > 10 4 K 
to cool at constant density for the duration of each timestep, as- 
suming a metallicity, Z = 0. Each gas particle with T < 2 x 10 4 K 
was then deemed eligible for star formation and we adopted the 
usual SFR prescription ^ = c»^p, where p* and p gas refer 
to the stellar and gas density respectively, f,/ v „ is the dynami- 
cal time of the gas and c, = 0.5 the star formation efficiency 
(this value is quite high to boost the SFR). Rather than spawning 
new (lighter) star particles, we implemented the above prescrip- 
tion in a probabilistic fashion. Assuming a constant dynamical 
time across the timestep, the fractional change in stellar density, 
Ap t /p* = 1 - exp(-cAf/?dyn)- We drew a random number, r, 
from the unit interval and converted a gas particle to a star if 
r < Ap*/p». 

After 1.5 hr x Gyr, a galaxy with stellar mass, M, = 4.5 x 
10 n M o formed with 75 per cent of the stars being born within 
the first Gyr. During the next 5 h~ l Gyr, we let the system re- 
lax to make sure that our new initial conditions are stable. It is 
worth mentioning that our aim here is not to perform a detailed 
study of the formation of an elliptical galaxy (which is poorly 
understood), but to construct a sensible starting point in order to 
assess the gross effect of forming the galaxy (i.e. from dissipa- 
tion) on our previous results, as well as the effect on the stellar 
distribution itself, which like the dark matter, is a collisionless 
fluid. 

From our new initial conditions, we perform the same simu- 
lation and analysis as described in the previous section. We use 
here A — 1 /2 r s and V gas = 320 km/s corresponding to the aver- 
age of the initial velocity dispersion of dark matter in the inner 
part. We justify shutting off cooling and star formation by assum- 
ing that the AGN activity balances cooling in our relaxed system, 
as suggested by Best (2006a). Although a significant fraction of 
gas has been converted into stars, the mass of the gas at 30h~ l 
kpc is still substantial (1O U M ) due to the subsequent inflow of 
hot gas from the outer part of the halo. 

Figure [5] shows the evolution of the DM density profile at 
t - 0,3,5, and 8/z~' Gyr. The DM cusp is again progressively 
flattened but a longer time is needed to obtain an effect of simi- 
lar magnitude to the model with no stars. However, no significant 
variation is found between 5 and 8 h Gyr suggesting that most 
of the effect happens during the first 5 h Gyr. It is also interest- 
ing to note that a similar effect is found with the stellar density 
profile, as shown in Fig. [6] In particular, we found a significant 
deficit of stars within 2 h~ l kpc. Here again, noticeable effects 
are found during the first 5 hr x Gyr. 
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Before drawing any conclusions, we need to check if other 
physical mechanisms are susceptible to play a role in the flatten- 
ing of the DM profile. In particular, after the phase of galaxy for- 
mation from baryonic collapse and/or sub-haloes merger, resid- 
ual clumps of gas are expected to exist and move through the 
virialized regions of the halo. If these structures are much more 
massive than DM particles and survive long enough from tidal 
interactions, then their orbital energy can be transfered to DM 
through dynamical friction. Consequently, the distribution of 
DM should be affected, especially in the inner part. This effect 
was clearly identified and characterized by El-Zant, Shlosman 
& Hoffman (TZOOlT) and El-Zant et al. d2004i l for both galaxy and 
clusters of galaxies. They found that the dynamical friction in- 
duced on clumps of gas with at least a mass of 0.01 percents of 
the total mass of the system can flatten the DM cusp after a few 
Gyr. 

In order to investigate if the mechanism proposed by El-Zant 
et al. (2001, 2004) is operating in our model, we have measured 
the dumpiness of gas present in our initial conditions. To do 
that, potential sub-structures have been identified by a friend-of- 
friend algorithm. We used different values of the linking length 
and, for each of them, we have estimated the number of sub- 
structures by using the following procedure. First, once a struc- 
ture is identified, the total energy of each particle is computed, 
with respect to the centre of mass, and those with positive energy 
are removed. The procedure is repeated with the new centre of 
mass, computed according to its usual definition, until no un- 
bound particles are found. Then, only structures with at least 57 
particles of gas, corresponding to objects with a mass greater 
than 10 times the DM particle mass (or 0.002 percents of the 
total mass of the system), were retained. 

We found that whatever the linking length, only one bound 
structure satisfying the above criteria is found at most (i.e 
the center of the halo itself). This proves that no bound sub- 
structures exist and that the gas is distributed in a smooth way 
(similar result is found for the stellar component). This is not 
surprising since we used a simplified prescription for the star 
formation in which the efficiency was quite high and no stan- 
dard criterion on the overdensity of the gas has been taken into 
account. Moreover, we let the whole system relax during 5 h 
Gyr after the merger event. Therefore, the lack of dumpiness 
of gas in our numerical model allows us to ignore the influence 
of dynamical friction. However, there is no denying that in re- 
ality, inhomogeneities and dumpiness of gas are expected to be 
present and so the dynamical friction acting on them may not be 
negligible. Assessing which mechanism will dominate is beyond 
the scope of the present paper as our aim was to provide a plausi- 
ble alternative mechanism, however we note that the effect found 
by El-Zant et al. (2004) at small fractions of the scale radius and 
after a period of 4 Gyr, seems to have a similar amplitude to our 
result. We can therefore reasonably think that some combination 
of dynamical friction acting on baryonic clumps and the mecha- 
nism proposed in this paper would significantly affect the slope 
of the inner DM profile. 

4. Discussion and Conclusions 

More and more observations indicate that the DM profiles of 
dwarf galaxies have a central DM core, in contradiction with 
results of numerical simulations. On the contrary, clusters tend 
to have a central cusp suggested by recent studies combining 
strong and weak lensing (see for instance, Limousin et al. 2006; 
Leonard et al. 2007; Umetsu, Takada, & Broadhurst 2007) even 
exception can be found (Sand, Treu & Ellis 120021 ). In this pa- 
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Fig. 5. DM density profile for V gas = 320 km/s and A = 1 /2 r s 
plotted at three different times, t — 3 (cyan), t — 5 (green) and 
t — 8 hT l Gyr (red). The black lines represent the initial DM pro- 
file (with a stellar component) which is steeper than our fiducial 
model discussed in section 2 and represented in red in the small 
panel. 
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Fig. 6. Stellar density profile for V gas - 320 km/s and A — 1 /2 r s 
plotted at t = (black), t = 3 (cyan), t = 5 (green) and t = 8 A -1 
Gyr (red).. 



per, we focus our attention on intermediate scales corresponding 
in particular to massive early-type galaxies. These are known to 
contain a central SMBH, whose mass scales with the spheroid 
velocity dispersion. The growth phase of this SMBH must have 
been nearly contemporaneous with the epoch of spheroid star 
formation. Which came first is controversial and not known, but 
it is eminently plausible that strong AGN activity occurs when 
the bulk of the spheroid stars have formed. For example, the 
epoch of activity of AGN as studied via hard x-rays yields a 
comoving growth rate that peaks at the same epoch as that of the 
cosmic star formation rate at z ~ 2. The most massive spheroids 
are most likely already in place. Nevertheless, the presence of 
SMBHs of mass ~ 10 9 M at high redshift (z ~ 6) presents some 
challenges to this theoretical model since it's not clear at all how 
can such massive objects and their host halos (with probably a 
mass greater than 1O 12 M ) can form so early in a ACDM cos- 
mology. One possible solution is that intermediate mass black 
holes form as remnants of Population III stars in the first gen- 
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eration of clouds (Spaans & Silk 2006 ). Then, accretion is the 
preferred mechanism given that the quasar luminosity function 
matches the local BH mass function for plausible accretion effi- 
ciency (Yu & Tremaine 2002 ). As far as the most massives halos 
being in place, yes, if they are rare as seems to be the case at high 
redshift. It's interesting to note that recent cosmological simula- 
tions from Li et al. 2007 of formation of high redshift quasar 
from hierarchical galaxy mergers seems to favor this scenario. 

Using N-body equilibrium systems, we have used a simple 
description of the AGN activity in which the energy transfer gen- 
erates a harmonic bulk motion on both sides of the gas core. We 
then investigate the impact of this mechanism on the DM core 
profile. For a typical giant galaxy of total mass 
we find that particularly strong effects occur when i) the initial 
velocity of the gas bulk motion is close to the velocity dispersion 
of the dark matter and ii) when the amplitude of the oscillations 
of the bulk motions extends over half the scale radius. In this 
case, the DM cusp progressively flattens after 4 - 5/i~' Gyr (with 
or without a stellar distribution) but cannot be well fitted by a 
Burkert profile. However, the effects obtained are not negligible 
and our model can be easily improved. For example, One could 
imagine that the direction and the efficiency of the outbursts do 
not remain constant in time. On the other hand, if we believe that 
ellipticals form by major mergers, then these merger events may 
lead to formation of black hole binaries. Binaries provide an al- 
ternative DM and stellar core heating mechanism. In this case, 
the combination of two AGN feedback modes may enhance the 
observed effects. Finally, the combination with other physical 
mechanisms may also improve the scenario. In particular, the 
dynamical friction acting on clumps of gas is expected to have a 
strong contribution, as shown by El-Zant, Shlosman & Hoffman 
(2001) and El-Zant et al. (120041 ). Taking into account this pro- 
cess, the effects obtained in the paper should be amplified. 

Our fiducial model is not unreasonable. For instance, Voit & 
Donahue (2005) claim that AGN outbursts of about 10 45 erg/s, 
lasting for at least 10 7 yr and occurring every 10 8 yr are required 
in order to explain the inner entropy profile in cooling-flow clus- 
ters. This equates to an energy outburst of 10 59 erg, comparable 
to the total kinetic energy periodically injected into the gas core, 
assuming V gns = 260 - 320 km/s. The frequency of the AGN 
emission is also in good agreement with the frequency v of the 
gas motions since v ~ 2.6 x 10 s /z 1 yr and v ~ 2.2 x 10 s h~ l 
yr respectively for V gas = 260 km/s and V gas = 320 km/s, us- 
ing A - 1/2 r s . Furthermore, it's worth mentioning that in the 
present case, the dynamical time of the core gas is comparable 
to the time between outbursts, taking a mean density of gas of 
~ 5 x 10 15 M su „.Mpc~ 3 in the central part. The displaced amount 
of gas is then expected to relax and sink toward the center be- 
tween 2 successive outbursts. Or, in a more realistic scenario, as 
gas is displaced (in a bi-conical nature), fresh gas can fall in per- 
pendicular to the jet axis on the same timescale. This can justify 
our assumption that the frequency of the AGN outbursts is equal 
to the frequency of the gas oscillations. Physically, this assump- 
tion seems reasonable since the two times are expected to be 
similar, as it is the accretion of core gas that drives the outbursts. 
However, we cannot totally exclude the possibility that either 
the length of an outburst is longer than the gas response time 
or the frequency of the outbursts is larger than the one used in 
our fiducial model. In both alternatives, the oscillations would be 
damped and the mechanism proposed inefficient. To finish, dark 
matter haloes and galaxies are continuously accreting matter. 
The corresponding gas merging/infall are responsible to drive 
the AGN activity which is proved to be sporadic with 10% duty 
cycle. Each outflow expels important amount of gas which be- 



comes a potential source of material for future BH accretion. 
According to these physical mechanisms, the periodical quasar 
activity is expected to be maintained during a long period, prob- 
ably over several Gyr, timescale required for our mechanism to 
operate. 

Our astrophysical justification for this timescale is the fol- 
lowing. Ellipticals most likely form by many minor mergers 
(Bournaud, Jog & Combes 2007). The major merger rate is too 
low. These are inevitably distributed over a few Gyr. If the gas 
supply in each merging event feeds the AGN as well as en- 
hances star formation, we inevitably expect the AGN activity 
and growth to be spread out, e.g. with a duty cycle of 10% over 
many dynamical times. The downsizing of AGN activity (and 
associated growth) supports such prolonged activity for typical 
ellipticals, with only the most massive forming early by major 
mergers. 

Moreover, it has been shown that the brightness profiles I{r) 
of nearly all ellipticals could be well fitted with Sersic / oc r 1/n of 
index n (see for instance, Caon, Capaccioli & D'Onofrio 1993). 
However, some recent studies have suggested that most giant el- 
lipticals exhibit "missing light" at small radii with respect to the 
inward extrapolation of the Sersic profile (see Graham [20041 or 
Ferrarese et al. 2006). Such observations are believed to be the 
consequences of dissipative galaxy formation mechanisms such 
as those involving AGN energy feedback. In particular, mas- 
sive SMBH binary tidal scouring to produce the observed mass 
deficit was advocated by Milosavljevic & Merritt (2003). We 
note here that AGN-induced gas bulk motions are equally likely, 
if not more natural, as a ubiquitous accompaniment to the gas- 
rich phase of massive spheroid formation. Our model assumes 
a substantial gas component not only during spheroid formation 
but also over the next Gyr of AGN activity. This constitutes a 
serious issue for all CDM models since gas infall is predicted. 
But as mentioned above, gas is expelled and or/induce to form 
stars at high efficiency as infrered from ultra-luminous infrared 
galaxies (ULIRGs) and submillimeter galaxies (SMGs). Wind 
are observed as well as high-efficiency star formation in these 
objects and an increasing frequency of AGN beyond z of 1 . It is 
at least 20% locally and given the different duty cycles for star 
formation and AGN. The feedback is likely to always be impor- 
tant. 

Indeed, in order to account for the correlation between 
SMBH mass and spheroid velocity dispersion, the period of 
SMBH growth by accretion (the case for accretion growth is 
powerfully argued by Yu & Tremaine 2002) must be followed by 
gas expulsion associated with SMBH-driven outflows. It is rea- 
sonable to argue that the expulsion phase is preceded by a phase 
when SMBH are subdominant with respect to the canonical re- 
lation but are still capable of driving bulk flows of the gas supply 
associated with sporadic, fueling-driven outbursts. Quenching of 
star formation by AGN activity, widely considered to play a role 
in massive elliptical formation, provides a plausible means of 
maintaining the gas bulk motions over several dynamical times. 
Since stars behave as a "collisionless" fluid similar to dark mat- 
ter, one should expect that AGN activities will affect not only the 
DM density profile but also the star density profile, confirmed by 
our results. In this case, the brightness in the inner part should 
be altered, in good agreement with observational studies such as 
Graham (2004) and Ferrarese et al. (12006b . 

Finally, we draw a parallel between giant galaxies and galaxy 
clusters. Indeed, Chandra observations have revealed the pres- 
ence of cold fronts in many galaxy clusters (see for instance 
Mark evitch et al. 200 0, 200 2; Vikhlinin, Markevitch & Murray 
1200 It Mazzotta et al. |2001| ). These cold fronts are supposed to 
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be due to the presence of a dense, cold gas cloud moving with 
respect to a hotter one. Some observations suggest that such cold 
clouds have "sloshing" motions around the center of the cluster, 
believed to be the consequences of either a recent major merger 
event or feedback processes from a central AGN. In fact, our 
model cannot be applied on such scales. Indeed, for a typical 
galaxy cluster of 10 15 h~ l M & , the velocity dispersion of DM in 
the central part is more than 1000 km/s. Then to obtain similar 
effects to our giant galaxy model, about 10 3 times more energy 
is required. It seems unrealistic that one or even a few AGNs lo- 
cated near the center of the cluster can provide such an amount of 
energy and induce gas bulk motions in a coherent way. However 
for massive galaxies, the situation is quite different. Not only are 
the energetics favourable, but the bipolar nature of the AGN out- 
burst must inevitably exert dynamical feedback during the gas- 
rich phase of spheroid formation. 
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